
The International Journal Of Engineering And Science (IJES) 
II Volume II 4 II Issue II 2 II Pages II PP. 19-3011 2015 II 
ISSN (e): 2319 - 1813 ISSN (p): 2319 - 1805 

Adaptive Grouping Quantum Inspired Shuffled Frog Leaping 

Algorithm 

12 2 2 

Xin Li , Xu Huangfu , Xuezhong Guan", Yu Tong , 

'School of Computer and Information Technology, Northeast Petroleum University, Daqing, China 
2 School of Electrical and Information Engineering, Northeast Petroleum University, Daqing, China 

ABSTRACT 

To enhance the optimization ability of classical shuffled frog leaping algorithm, a quantum inspired shuffled 
frog leaping algorithm with adaptive grouping is proposed. In this work, the frog swarms are adaptive grouped 
according to the average value of the objective function of child frog swarms, the frogs are encoded by 
probability amplitudes of Multi-Qubits system. The rotation angles of Multi-Qubits are determined based on the 
local optimum frog and the global optimal frog, and the Multi-Qubits rotation gates are employed to update the 
worst frog in child frog swarms. The experimental results of some benchmark functions optimization shows that, 
although its single step iteration consumes a long time, the optimization ability of the proposed method is 
significantly higher than the classical leaping frog algorithm. 
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I. INTRODUCTION 

Shuffled frog leaping algorithm (SFLA) is a new heuristic cooperative search algorithm, which 
simulates the foraging behavior of a group of frogs jumping in wetlands [1]. As a new bionic intelligent 
optimization algorithm, SFLA combines the advantages of memetic evolutionary algorithm and particle swarm 
optimization algorithm, with the concept of simple, less adjustment parameters, calculation speed, strong global 
search capability, easy implementation, currently used in the rough set attribute reduction [2], fuzzy-means 
clustering [3], drawing water bath to control [4], cooperative spectrum sensing cognitive radio [5- 6], speech 
recognition [7], etc. In terms of performance improvements SFLA, paper [8] by introducing in the frog 
population attraction and repulsion mechanism, effectively avoiding premature convergence; paper [9] by 
introducing the extreme dynamic optimization to improve the optimization efficiency of the algorithm; paper 
[10] through the integration, simulated annealing, immunizations, Gaussian mutation, chaotic disturbance, 
enhanced algorithm optimization capabilities; paper [11] by using a dispersion and fitness, improve the 
convergence speed and accuracy. However, these improvements do not consider the effect of the grouping 
SFLA for optimal performance. Since SFLA each group to update only a worst frog, in general, when a fixed 
number of frog, the fewer packets (the more the number of the group frog), the higher the computational 
efficiency, and optimization capability is relatively weak, so frog group grouping has an important influence on 
the optimization of the performance of SFLA. Taking into account the optimal number of packets SFLA usually 
related to specific issues as well as the optimization phase, this paper presents an Adaptive Grouping Quantum 
Inspired Shuffled Frog Leaping Algorithm. Quantum computing is an emerging interdisciplinary; combining 
information science and quantum mechanics, and its integration with intelligent optimization algorithms began 
in the 1990s. Currently, there are a lot more mature algorithms, such as quantum-behaved particle swarm 
optimization algorithm [12], quantum inspired evolutionary algorithm [13], quantum derivative harmony search 
algorithm [14], quantum inspired immune algorithm [15], quantum inspired genetic algorithm [16], quantum 
inspired differential evolution algorithm [17]. In addition to the literature [12] using the real-coded, the rest are 
encoded using a single bit probability amplitude. Single-bit probability amplitude coding disadvantage is that 
the adjustment of a qubit can only change one gene locus on the individual, while the probability amplitude for 
multi-bit coding, adjust a qubit, can change all of the ground state probability amplitude in multi-bit quantum 
superposition states , thereby changing the position of all genes on the individual. Therefore, this paper will also 
propose a new multi-bit probability amplitude coding mechanism to further improve the SFLA of optimizing 
capacity and optimize efficiency. The results of the optimization criterion function demonstrate the superiority 
of the proposed method. 
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II. SHUFFLED FROG LEAPING ALGORITHM 

First, There is N frogs x > ~ *- x n' x t*' ' x t«' , calculated the target value f^ x t^' ' 2 ' '^of each 
frog, sort all Target value from good to bad. Order N = k x m, k is the number of sub-group, m is the number of 
subgroups frogs, the N sorted frogs cycle divided into k sub-groups, namely: the first sub-group divided into 
optimal frog, The second sub-group divided into Suboptimal frog, and so on, until all the frogs have been 

allocated. In each subgroup, the best and worst frogs were recorded as * and » , the whole frog group of the 

X x 

best frog recorded as * , in each iteration, the following formula is updated sub-group » . 

D = r K - w „) (1) 

< = *» + " (2) 

where ° m " ~ D > ~ D "" , ^ ' is the /* value of the vector D , i = 1, 2,..., n, D '■»" is the biggest leap frog step, 
r is a random number between 0 and 1 . 

If x - is better than x » , then x » instead of"*" ; otherwise, Xb instead of s , repeat the Eq.(l) and (2). If the x * 

inferior to" 1 " ,, then randomly generate a solution to replace X " . So the cycle until the termination condition is 
met. 

III. ADAPTIVE GROUPING STRATEGY 

In SFLA, by updating the group worst frog to achieve local optimization, by mixing and re -grouped in 
each group to achieve global optimization. From the optimization mechanism, the group number reflects the 
optimization of the global characteristics, while the number of frogs in the group reflects the optimization of the 
locality. So, how to achieve the balance of global and localized search, is worthy of in-depth study of the 
problem, and this balance by grouping decision. However, the balance of global and localized search is often 
also associated with specific issues and optimization phase, grouping of frog group cannot be fixed pattern. 
Therefore, this paper proposes an adaptive grouping strategy. 
Let the total number of frogs be N, we rewritten N as following equation 

N = n l x = n 2 x m 2 = — n y. m ^ 

where 12 s 1 2 

Therefore, there are s kinds of groupings. Namely: group, each group has m ' frogs; " 2 group, each group 
has ™ 2 frogs;...; « group, each group has • frogs. As an example of the minimum value optimization, 

/ n 

recorded the average value of the objective function as ""' , Frog group is divided into 1 groups, each with 
frogs, the average value of the * group objective function is " s , the " s number is "* , the 

/' > / n 1 

" g number is " . Adaptive grouping strategy proposed in this paper can be described as follows. 

^ «. \n, t . => n. \n < nl \n,. , n, 



, k < s , I m , , => m , T - 1 k > 1 , m, , => m, 

(1) If 1 ,than 1 i + 1 * , (2) If L , than 1 4 . 

For this strategy, we explain below. n ^> n " represents the sub-groups of "* less than have advantages, 
This means that the optimization approach tends more parallel sub-groups, and therefore need to increase the 

number of sub-groups; Conversely, < " t2 represents the sub-groups of less than " K have 

disadvantages, this means that the optimization approach tends to parallel the minority carrier group, hence the 
need to reduce the number of sub-groups. 

IV. MULTI-BIT QUANTUM SYSTEM AND THE MULTI-BIT QUANTUM ROTATION GATE 
4.1. Qubits and single qubit rotation gate 

What is a qubit? Just as a classical bit has a state-either 0 or 1- a qubit also has a state. Two possible 

states for a qubit are the state ' °^ and ' , which as you might guess correspond to the states 0 and 1 for a 
classical bit. 

Notation like 1 * is called the Dime notation, and we will see it often in the following paragraphs, as it is the 
standard notation for states in quantum mechanics. The difference between bits and qubits is that a qubit can be 
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in a state other than l0> or ' . It is also possible to form linear combinations of states, often called 
superposition. 

I (j>) = cos 6 I 0) + sin 9 1 1) = [cos 8 sin6>] T , 

where 9 is the phase of , cos e and Sln 9 denote the probability amplitude of . 

In the quantum computation, the logic function can be realized by applying a series of unitary transform to the 
qubit states, which the effect of the unitary transform is equal to that of the logic gate. Therefore, the quantum 
services with the logic transformations in a certain interval are called the quantum gates, which are the basis of 
performing the quantum computation. A single qubit rotation gate can be defined as 

[cos A # - sin A 6> ] 
R(A0) = \ 

I sin A 0 cos A 8 J 

fcosfll |~cos(0 + A0)l 

1^ = 1 i R(A0) = \ 

sin 0 1 ,\d>) , c , . sin(0 + AO) T . 

Let the quantum state L J , and r can be transformed by L J . It is obvious 

that R(A0) shifts the phase of W 

4.2. The tensor product of matrix 

Let the matrix A has m low and n column, and the matrix B has p low and q column. The tensor product 
of A and B is defined as. 



A ® B = 



A . 

where ''' is the element of matrix A. 







L 


A ,„ B 1 




A 22 B 


L 


A, B 


1 M 


M 


M 


M 


1 


A m , B 


L 


A B J 



(6) 



4.3. Multi-bit quantum system and the multi-bit quantum rotation gate 

In general, for an n-qubits system, there are 2" of the form ' x,XlL x ^ ground states, similar to the 
single-qubit system, n-qubits system can also be in the a linear superposition state of 2" ground states, namely 

i i i 

'A^L $U=Z Z L Z x. ljC l*3 L *»> = KoLO fl 00Ll L "llLl^' 

where *'* 2L *« is called probability amplitude of the ground state 12 » , and to meet the following 
equation. 

z z l ik, „r=i- 

*> (8) 
Let _ cos # 1°) + sin .1 > ^ accor( j m g to the principles of quantum computing, the 1 ^1^= L r»> can be written 



[cos £? ] f cos # 1 

> =1 </> )® I 0,) ® L I > = I ® L ® I 

L sin °» J 

I cin # cin f) T cin # I 

(9) 

It is clear from the above equation that, in an n-qubits system, any one of the ground state probability amplitude 

(0 0 L 0 \ 0 
is a function o/n-qubits phase 1 ' 2 ' ' » ' , in other words, the adjustment of any ' can update all 2" 

probability amplitudes. 

In our works, the n-qubits rotation gate is employed to update the probability amplitudes. According to the 

principles of quantum computing, the tensor product of n single-qubit rotation gate ' is n-qubits rotation 
gate. Namely 
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R(A0.A0 2 L A0J = «(A0 ( ) ® «(A6> 2 ) ® L ® R(A 0 : ), 

[cos A 6> - sin A 0. ~| 
fi(A9.)= ' I, / = 1,2,L ,n. 

, sin AS cos A ^ 

where L ' ' J 

Taking n=2 as an example, the R< - A0 i AS 2^ can b e rewritten as follows. 

[cos A 0 [ cos A 9 2 - cos A 0^ sin A 9 2 - sin cos A9, sin A # t sin A $ 2 
J cos A 0 l sin A 0, cosA^cosAS, - sin A 0 t sin A sinA^cosAi?, J 

lsinA#cosA# -sinA#sinA# cosAff cosAfl - cos A 0 sin A 0 I 

I I 
[ sin A 0 { sin A # 2 sinA#,cosA# 2 cos sin Afl, cosA^cosAff, J 

It is clear that 

R ii (A0 1 A0 2 h A0 n )\<t> [ 0 2L 0J=\j l )®\0,)®L ®\<pj, 
where = cos( ^ + Aff ' ) 1 0) + sin( ^ + A6) /> u >- 



(10) 



R(A( 



(11) 
(12) 



V. HUFFLED FROG LEAPING ALGORITHM ENCODING METHOD BASED ON MULTI- 
BITS PROBABILITY AMPLITUDES 

In this paper, the frogs group is encoded by multi-qubits probability amplitudes. Let N denote the 
number of particles, D denote the dimension of optimization space. Multi-qubits probability amplitudes 
encoding method can be described as follows. 

5.1. The number of qubits needed to code 

For an n-bits quantum system, there are 2n probability amplitudes, which can be used directly as a 

result of an individual encoding. In the D-dimensional optimization space, it is clear that D - 2 . Due to the 

constraint relation between each probability amplitude (see to Eq.(10)), hence D < 2 . For the D-dimensional 
optimization problem, the required number of qubits can be calculated as follows. 

n = log(D) + 1 . 

5.2. The encoding method based on multi-qubits probability amplitudes 

First, generating randomly N n-dimensional phase vector S ' '■ , 1 ~ 1 ' 2 ' L ' N ' as follows 



(13) 



(14) 



" :" 

0.. = In x rand 7 = 1 2L n 

where '' , ran d is a random number uniformly distributed within the (0,1), 

I </>. > = cos0.. I 0) + sin 0 . 1 1) 

Let u " , Using Equation (11), we can obtain following N n-qubits systems 

''»^ L , '^,^2 L >->>j !^Ai L . i n each of the quantum system, the first D probability 
amplitudes can be regarded as a D-dimensional particle code. 

VI. THE UPDATE METHOD BASED ON MULTI-QUBITS PROBABILITY AMPLITUDES 

In this paper, the multi-bit quantum rotation gates are employed to update particles. Let the phase 

F = [9 .,0 ,,L 0 ] 

vector of the global optimal frog be 8 s 8 s " , the phase vector of the group optimal frog is 
F b " [^6i»^62' L ' S i„J j fjjg phase vector of the group worst frog F " ~ ^»i'^»» ,L ■^»»1 . Similar to traditional 
SFLA, for each subgroup, we only just need to update 

By formula (9) it is clear that, once » has been updated, all its corresponding probability amplitudes will be 

F 

updated. To improve the search capability, in an iteration, all phases - are updated in turn, which allows all 

A 9 

particles are updated n times. Let 0 denote the phase update step size, the specific update can be described as 
follows. 
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Stepl. Set y=l, " L " 12 " " ' 2 "' . 

Step2.Set A ^= A ^= L = A ^=°. 

Step3: Determine the value of the rotation angle, where the sgn donates the symbolic function. 

\0\. - 0.. I < it A 9 b . = sgn(0* - 0. )A8 a 

If "■' '■' , then ' J 0 " 0 . 

10,'. - (9.1 > >r A6>* = - sgn((9' - 0„)A0 n 

If '" " , then ,J " J 8 0 . 

l0,-0„l>>r Afl* = sgn(<9 , - 6> )Atf„ 

If 8j ' J , then " " ,J 0 . 

\9 . - <5U > it A0 S = -sgn(0 . - 0 )A8 n 

If SJ " ,then ' J SJ ' J °. 

Step4 A9 ' = °' 5A9 "i + °- 5A ^J = RJA0 l ,A0 z ,L ,&0JF w (,0) 

c , c t A9 = rnds ' A9 Q , . F (2) (8) = R (A 8, , A 0 7 ,L , A 0 )F "' (0) T „ 

Step5: Let J , r " fl 5 is a random number of -1 to 1. " " 1 2 " » v .If 

F'"(0) ■ u „ F W (0) F (0) = F O) (0) ■ F (0) = F ,2] (0) 

is better than » , then 11 » , otherwise, » » v . 

Step6: If < " , thenj = j + 1, back to step2, Otherwise, end. 

VII. ADAPTIVE GROUPING QUANTUM INSPIRED SHUFFLED FROG LEAPING 

ALGORITHM 

Suppose that, Af denote the number of frogs, D denote the number of optimization space dimension. 

11 » » , ' and ' is positive integers and 12 » 1 2 * . For adaptive 

grouping quantum-inspired shuffled frog leaping algorithm, called AGQISFLA, the optimization process can be 
described as follows. 
(1) Initialize the frog group 

According to Eq.(12) to determine the number of qubits n, according to Eq.(13) initialize phase of each 
particle, according to Eq.(9) to calculate the probability amplitude of 2" each particle, where the first D 

probability amplitudes are the coding of the particles. Set the / h probability amplitude of the z* particle be '' , 
coding result can be expressed as the following equation. 



P, = [x 2l ,x 22 ,h ,x 2D ] 
L L L L L L L 

k =[*»,. *,a. L < X »J 



(15) 

A 0 

Initialization phase update step 0 , the limited number of iteration G. Set the current iteration step f=l. 
(2) Calculation of the objective function value 

[MinX ,MaxX .] * 

Set the y'-dimensional variable range be ' ' , because of the probability amplitude " values in 

the interval [0,1], it is need to make the solution space transformation. The transformation equation is below. 



1 

X =—[MaxX (l + x ) + MinX (1 - x )], 

2 (16) 

where i= 1, 2, —,N,j = 1, 2,—, D. 

With the above formula, calculate the objective function values of all frogs. And ascending objective function 

F = [0 .,0 ,,L ,0 ] f 
value, global optimal frog phase be s s s s " , global optimal objective function value be ! ,tLS is 

the total number of iteration. 

(3) Frog Segmentation 

The sorted frogs cycle into L sub-groups, each group have M frogs, namely: Optimal frog divided into the 
first subgroup, suboptimal frog divided into the second subgroups, and so on, until all frog allocated. 

(4) Update subgroup worst frog 

The Each subgroup evolution times is LS, in each evolution, update subgroup worst frog , follow the steps in 
the previous section (1) to step (6). As such, the worst frog of each subgroup are updated LS x n times. 

(5) Adaptive calculation of the number of subgroups 
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/' 

In each sub-group, the average value of the objective function is " n , where /= 1, 2,.., L. Mixing the 
subgroups, and arranged in ascending order according to the value of the objective function, remember the 

/ f' < f n 

objective function is the average of the entire frog groups is " " , the number of °' s " s is 1 , the number 

f' > f n 1 

of " g is L .The following two steps to recalculate the number of sub-groups are L and the number of 

subgroups frog are M. 

n\ > n\ [L = l M \n[ < n L [l = l k _ { 

(a)If L , than 1 t+l . (b) If 1 ,than 1 '-' . 

(6) Update the global optimal solution 

Let the optimal frog phase be 8 81 ' 82 ' ' 8 " , the corresponding objective function value be 8 . If 

f < f f = f F = F f = f F = F 

1 1 , then 1 * , s 8 ; otherwise " 8 , * 8 . 

(7) Determine the termination condition 

If t <G, t = t + 1, back to (3); otherwise, save the result, end. 

VIII. COMPARATIVE EXPERIMENT 

In this study, the 25 standard test functions are employed to verify the optimization ability of 
AGQISFLA, and compare with the traditional leapfrog algorithm (SFLA), adaptive grouping leapfrog algorithm 
(AGSFLA). All test functions are minimal value optimization, All functions belong to minimum optimization, 

where D is the number of independent variables, Q is the solution space, x is the exact minimum point, 

f ( x * ) 

is the corresponding minimum. 
8.1. Test function 

(1) /, (X ) = £ ° * ( S ; Q = [-100,100]°; X ' = [0, 0,L , 0]; f(x') = 0. 

(2) f 2 (X) = £ °\x. \ + Yl °, k p fl = [-100,100]°; X " = [0,0, L ,0]; f(X*) = 0. 

(3) = £ D ,(X '-i*p 2; = [-100,100]° ; X " = [0,0, L ,0]; f(x") = 0. 

(4) f(X ) = max(|.v I); Q = [-1 00,1 00]° ; X " = [0,0, L , 0]; f(x') = 0. 

l<i<D 1 ' 1 

(5) f 5 (X ) = £ D ^(100(x : + I - x. 2 ) 2 + (x. - l) 2 ); Q = [-100,100]"; X " = [0, 0,L , 0] ; f(x") = 0. 

(6) f 6 (X ) = £ " [x, + 0.5] 2 ;H = [-100,100]°; X + = [1,1, L ,1]; f(X + ) = 0. 

(7) / 7 (X) = °_ 1 ix *( l + random (0,1)); fi = [-100,100]°; X " = [0, 0,L , 0]; f(X*) = 0. 

(8) f(X ) = [x 2 - 10cos(2^-jc) + 10]; fl = [- 1 00, 1 00]° ; x'=[0,0,L ,0]; /(X*) = 0. 



(9) / 9 (Z ) = -20 exp(-0.2^|— ^ " ^x 1 : ) - exp( — £ ° | co s( In x. )) + 20 + e; 

Q = [-100,100]°; X* = [0,0,L ,0]; f(X") = 0. 

(10) f w (x ) = — ^ Y " x; -rr ° cos(^) + i ; a = [-100,100]°; jt* = [o.o.l ,0]; /(x*) = o. 

4000 '■' '"' v,- 

(11) f u (X) = —V ° (-v 4 -16j; 2 +5A:.) + 78.3 323314;fi=l-100,100] D ;A:* = -2.903534;/(A'*) = 0. 

D 

(12) /.2< x > = -^-L10sin 2 (jr y,) + £ - D'U + 10 »in ! (jr + <?„ - D* ] + £ *,«<*,, 1 0, 1 00.4 ); 

f & ( x — a~)'" , x > <a 

' i 

u{x ,a,k,m) = < O, — a <. x a y — 1+ — ( x + 1 ) ; 

i » ' ' 4 ' 

^ /: ( — x . — (7 ) , x t , — a 
a = [-100,100]°;^ ' = [-1,-1, L ,-l];/(X*) = 0. 
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(13) /,,(X ) = — [sin 3 (3tf*,)+ V ""'(*, - 1) 2 (1 + sin 2 (2;r*. +| )) + 

D '"' 

(x D - 1) 2 (1 + sin 2 (2^-.Y D ))] + Yj ° «(*,.5,100,4); 

I /: ( x . — a )" , x > a 
I 

M(x.,CE,«,m) = -j 0, — a < x . < a 
\k(—x j — a )'" , x. , — a 

Q = [-100,1 00] D ;X* = [-1,-1, L ,-l];/(X + ) = 0. 

(14) f l4 (X ) = £ °~ (** + 2^. 2 +| - 0.3cos(3^x.)cos(4^-x. + | ) + 0.3); 

= [-100,100]°; X* = [0,0,L ,0]; /(X*) = 0. 

2 

(15) /„(*) = Z° ,Z J- cos(y is ) + l); 

" '-' 4000 

=ioo(*. - x)y + d- x/- 

SI = [-100,1 00] D ;X * = [l.l.L ,1]; /(X*) = 0. 

(16) / 16 (X)= 2 + 10 - l ' Ji -2) 2 + 5 (* 4i -, - ^4,) 2 + 0 4i - 2 - 2 ^4,-,) 4 +1 °( JC 4,-3 + 10X 4| ) 4 ]; 

fl = [-100,100]°; X*=[0,0,L ,0]; /(X*) = 0. 

(17) /„(X)= ^ "j lg(jt < ,JC ^ )+ S(* D >*i)s(- t '>') = (** + y 2 )°' 25 x[sin 2 (50(x 2 + /) 0 '') + l]; 

= [-100,100]°; X* = [0,0, L ,0]; /(X*) = 0. 

k | < 1 / 2 

[round (2x.) 12, > 1 / 2 

i2 = [-100,100]"; X* = [0,0,L ,0]; /(X*) = 0. 

(19) / 19 (X ) = £ *"7 [a " cos(2«-i'U f + 0.5))]} '""(a* cos (^ fot )); 

a = 0.5;fo = 0.3; ft max = 30; Q = [-100,100]° ; X* = [0,0, L 0]; /(X*) = 0. 

D D D 

(20) / 20 (X ) = ^ a- 2 + 0.5/jc.) 2 + 0.5«.) 4 ; = [-100.100]°; X * = [0, 0,L 0]; /(X + ) = 0. 



(18) / lg (X ) = 10D + £ , t (y' - 10cos(2^y.)); 



(20) 




D 

= I *, 


(21) 


/ 2I (X) 


= !,. 


(22) 


f 22 (X) 


D 


(23) 




Z) -1 

= - 1 






Si = [ 



-(* +x +0.5xx ) fl ; 

exp ) x cos(4a/jt : . + x. +1 + 0.5x.x f+1 )J + D 



(24) / 24 (X ) = 1 - exp(-0.5^ " * 2 );fi = [-1,1]° ; X + = [0,0, L ,0];/(X + ) = 0. 

(25) / 2 ,(X ) = 1 - cos(2^-^^ ^ a: 2 ) + 0.1^ ^x 2 ;^ = [- 1 00, 1 00] ° ; X + = [0,0, L ,0];/(X + ) = 0. 

8.2. The experimental scheme and parameter design 

The total number of frog group is set to N = 100, in order to achieve adaptive grouping, the N is 
decomposed into the following eight situations, as follows. 

iV = 2 x 50 = 4x 25 = 5 x 20 = 10x 10 = 20x 5 = 25 x 4 = 50 x 2 = 100 x 1 ryj\ 

Among them, each instance is a grouping scheme, the first number is the number of sub-groups, and the second 
number is the number of frogs in the group. 

SFLA and AGSFLA limit the number of iteration steps is set to G = 100 and G = 1000, Biggest jump step taken 
D »" ~ 5 ; limited iteration number of AGQISFLA is set to G = 100, The subgroups iterations of these three 
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algorithms is set to L5'=15;phase update step is set to AS ° ~ °- 057r 

For SFLA, each function are used eight kinds of groupings were optimized 50 times, namely, each function is 
independent optimization .400 times. 400 times optimal results averaged and the average value of a single 
iteration of the running time as a comparative index; for AGSFLA and AGQISFLA, The initial number of frog 
groups is L = 10, the group number of frog is M = 10, Each function independently optimized 50 times, taking 
the average of 50 times optimal results, with the average value of a single iteration of the running time as a 
comparison index. 

8.3. Comparative Experiment Results 

Experiments conducted using Matlab R2009a. For comparison, the average time of a single iteration of 

f T 0 

the function J ' is set to > , the average optimal results for ' , i = 1, 2,... ,25. All test functions, Taking G=50 as 

f 

an example ( 16 for D = 52), the results of such comparison are shown in Table 1, he average optimization 
results for D=100, are shown in Table 2. 

f T A T F 

For the function J - ', the average time of the single iteration of SFLA, AGSFLA, AGQISFLA is set to 1 , '' , 

T a 0 A 0 F 0° 

' .the average optimal results is set to : , ' , ' .To facilitate a further comparison, Taking AGSFLA and 
SFLA as example, there are the following two formulas. Eq.(18) is the ratio of the average running time. Eq.(19) 
is the ratio of the average optimal results 

T F ' 25 (18) 

a" _ g "e," ' Q, F 

Q' ' 25 (19) 

For these three algorithms, the ratio of average running time and the average optimal results is shown in Table 3. 
Table 1 The 50 times contrast of the average results for three algorithms optimization (D = 50) 



SFLA 


AGSFLA 


AGQISFLA 




r,(j) 


0, 


r,(i) 


o, 


r,(») 


O, 


G = 1 0 3 


G = 1()' 


G =10' 


G = 10 3 


G = 1(T 


/, 


0.0146 


2.03E+03 


3.09E+02 


0.1452 


2.61E+02 


2.59E-008 


0.1700 


3:26E-015 


/, 


0.0163 


4.23E+08 


3.41E+02 


0.1103 


5.30E+05 


5.64E+02 


0.2091 


5:19E-008 


f, 


0.0984 


5.16E+04 


8.11E+03 


0.5064 


5.55E+03 


3.35E+02 


1.2738 


L17E-015 


ft 


0.0190 


14.9593 


10.4918 


0.1055 


15.3773 


12.8575 


0.2239 


0:514605 


f, 


0.0372 


7.04E+07 


7.73E+05 


0.2028 


7.42E+05 


2.03E+02 


0.3585 


48:24962 


ft 


0.0279 


1.58E+03 


43.23333 


0.1434 


85.70000 


34.30000 


0.2793 


0 


fi 


0.0341 


2.84E+07 


3.24E+03 


0.1532 


4.92E+04 


0.001770 


0.2827 


3:55E-028 


f, 


0.0249 


2.40E+03 


9.78E+02 


0.1377 


1.09E+03 


9.71E+02 


0.2221 


3:16E-013 


f, 


0.0316 


17.18656 


16.07727 


0.1465 


14.31819 


13.61492 


0.4036 


0:1850413 


f 

J 1 0 


0.0308 


1.62E+02 


3.149318 


0.0740 


4.434261 


0.002976 


0.3502 


L16E-015 


fu 


0.0317 


1.35E+04 


1.21E+02 


0.1606 


1.45E+02 


12.72304 


0.3003 


6:1248557 


fz 


0.0773 


8.40E+06 


11.55566 


0.4046 


2.94E+02 


17.54648 


0.6654 


2:99E-005 


As 


0.0772 


2.59E+07 


8.99E+03 


0.3996 


5.32E+04 


1.00E+02 


0.7448 


0:0322664 


fu 


0.0366 


6.05E+03 


9.57E+02 


0.1735 


7.94E+02 


11.75275 


0.4062 


L42E-013 


fx, 


0.3756 


3.33E+13 


2.35E+09 


0.7196 


1.41E+10 


2.52E+03 


3.1635 


2:90E+002 


fe 


0.0515 


3.76E+07 


3.98E+05 


0.2787 


1.97E+05 


2.26E+03 


0.6207 


9:26E-023 


/,, 


0.0515 


1.93E+02 


1.62E+02 


0.3485 


1.77E+02 


1.65E+02 


0.5558 


0:0016304 


f, 


0.0447 


2.42E+03 


9.61E+02 


0.3161 


1.24E+03 


9.68E+02 


0.3490 


L26E-012 


fit 


0.4529 


58.03894 


41.94107 


0.8332 


60.97842 


44.34244 


3.7331 


6:8459530 
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Ao 


0.0202 


1.55E+07 


1.19E+04 


0.1772 


3.08E+04 


/.ooh+03 


0.2160 


C.Jlc A1C 

5:43E-U15 


A. 


0.0338 


4.17E+65 


2.50E+47 


0.2809 


4.13E+41 


1.53E+15 


0.3652 


l:27E-206 


A* 


0.0167 


1.47E+02 


74.83271 


0.0969 


71.211607 


46.331449 


0.1750 


6:76E-009 


/,, 


0.0348 


44.76880 


43.062524 


0.1367 


45.373445 


42.675330 


0.2528 


0:1789964 


/« 


0.0158 


0.073135 


0.0317157 


0.0996 


0.0041286 


1.01E-14 


0.1133 


0 


A, 


0.0173 


1.92E+02 


15.005397 


0.1037 


24.769981 


4.357839 


0.1543 


4:72E-014 


Table2 The 50 times contrast of the avera 


ge results for three algorithms optimization (D = 100) 


/, 


SFLA 


AGSFLA 


AGQISFLA 


r (.s) 




T, ( s ) 




T. ( j ) 


0 : 


G = 1 0 2 


G = 1 0 3 


G = 1 0 2 


G = 10' 


G = l o 3 


/, 


0.0226 


4.99E+03 


1.02E+03 


0.2278 


2.16E+03 


0.061430 


0.2606 


8:91E-015 


A 


0.0246 


3.64E+51 


6.45E+04 


0.1723 


1.27E+35 


1.15E+03 


0.3227 


l:15E-007 


A 


0.1544 


2.10E+05 


3.12E+04 


0.7998 


1.41E+04 


5.87E+03 


1.9394 


1:61E-013 


A 


0.0292 


18.15317 


13.80720 


0.1683 


17.63319 


11.72408 


0.3488 


0:559791 


A 


0.0589 


1.48E+08 


2.98E+06 


0.3147 


2.58E+07 


4.78E+02 


0.5650 


69:8203 


f. 


0.0442 


3.68E+03 


3.46E+02 


0.2277 


8.68E+02 


3.80E+02 


0.4337 


0 


A 


0.0531 


1.43E+08 


2.38E+05 


0.2324 


6.12E+06 


6.29E+03 


0.4421 


l:39E-027 


A 


0.0385 


5.66E+03 


2.59E+03 


0.2200 


4.01E+03 


2.95E+03 


0.3377 


013 


A 


0.0502 


18.44872 


17.75022 


0.2237 


17.56241 


12.90171 


0.6207 


0:404291 


/,. 


0.0487 


3.37E+02 


6.909176 


0.1128 


27.02232 


0.077702 


0.5527 


9:11E-014 


/,, 


0.0495 


1.62E+04 


2.06E+02 


0.2549 


5.62E+02 


13.19321 


0.4708 


9:089264 


f» 


0.1222 


1.80E+07 


17.88713 


0.6367 


3.51E+05 


16.03746 


1.0287 


l:23E-004 


/,3 


0.1208 


5.76E+07 


7.37E+04 


0.6048 


1.03E+07 


2.06E+02 


1.1594 


5:01E-002 


/» 


0.0561 


1.48E+04 


3.04E+03 


0.2604 


5.33E+03 


32.15902 


0.6415 


2:04E-013 


/is 


0.5742 


1.95E+14 


2.05E+10 


1.1437 


6.90E+12 


1.47E+04 


4.7639 


3:99E+002 


/,« 


0.0790 


5.23E+07 


7.03E+05 


0.4235 


1.03E+06 


1.72E+04 


0.9684 


4:32E-022 


/l, 


0.0800 


4.16E+02 


3.54E+02 


0.5331 


3.73E+02 


1.43E+02 


0.8364 


0:091603 


/,. 


0.0703 


5.37E+03 


2.70E+03 


0.4950 


4.30E+03 


1.88E+03 


0.5379 


5:11E-012 


/,, 


0.6933 


1.39E+02 


1.16E+02 


1.2734 


1.44E+02 


21.29647 


5.7135 


9:935560 


A, 


0.0319 


1.17E+05 


4.09E+04 


0.2741 


1.15E+05 


1.89E+04 


0.3428 


3:88E-014 


/,, 


0.0526 


7.74E+84 


1.71E+55 


0.4231 


1.87E+60 


1.63E+17 


0.5483 


9:36E-187 


A J 


0.0256 


3.56E+02 


1.98E+02 


0.1549 


2.64E+02 


1.58E+02 


0.2759 


8:34E-008 


As 


0.0546 


93.59604 


92.42014 


0.2130 


94.61277 


90.05981 


0.4037 


0:9899643 


A 4 


0.0245 


0.183484 


0.096426 


0.1536 


0.029085 


1.12E-08 


0.1811 


0 


A, 


0.0267 


4.92E+02 


91.40583 


0.1608 


2.06E+02 


14.42596 


0.2436 


2:67E-013 



Table 3 the ratio of average running time and the average optimal results 



D 


T IT r 


Qt i Q* 


T o /T , 




T Q IT 1 




G = 10 2 


G = 10 3 


G = 10 2 


G = 10 3 


G = 10 2 


G = 10 3 


50 


5.495 


0.2784 


0.5133 


2.107 


0.0819 


0.0418 


10.06 


0.0067 


0.0111 


100 


5.486 


0.4143 


0.3492 


2.117 


0.0602 


0.0567 


10.04 


0.0054 


0.0081 


AVG 


5.491 


0.3463 


0.4312 


2.112 


0.0710 


0.0492 


10.05 


0.0060 


0.0096 



From Tab. 1 -Tab. 3, the introduction of adaptive grouping strategy and quantum computing makes AGQISFLA 
single-step running time of about 10 times that of traditional SFLA. Therefore, in order to enhance the fairness 
of the comparison results, we not only need to look at the same contrast iteration number, and must be further 
investigated algorithm to optimize the results of comparison in the same time. This is the fundamental reason for 
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the SFLA and AGSFLA the iteration number is set to G = 100 and G = 1000. According to the optimization 
results of fl ~ f25, when G = 100, AGSFLA is about 0.35 times of SFLA; when G = 1000, AGSFLA is about 
0.43 times of SFLA. This shows that the introduction of adaptive grouping strategy can indeed enhance the 
ability to optimize the algorithm. When G = 100, AGQISFLA (G = 100) is about 0.07 times of AGSFLA; when 
G = 1000, AGQISFLA (G = 100) of approximately 0.05 times of AGSFLA. This shows that the use of multi-bit 
probability amplitude coding and evolutionary mechanisms can indeed improve the algorithm optimization 
capabilities. From 



Table 3, in the same iteration steps, optimization results of AGQISFLA only 6/1000 of SFLA; at the same time 
optimization, optimization results of AGQISFLA only one percent of SFLA. Experimental results show that the 
adaptive grouping and multi-bit probability amplitude coding can indeed significantly improve the ability to 
optimize the traditional leapfrog algorithm. 

In this paper, when D = 50 and D = 100, the average results of the algorithm AGQISFLA is show in Fig.l, the 
average results of 50 times each function after optimization, the figure shows that when the dimension of 
AGQISFLA increases, there is a good stability. 




Figure 1 When D = 50 and D = 100, contrast optimization results of AGQISFLA 

8.4 Analysis of experimental results 

About adaptive grouping strategy, when fl ~ f25 is 100-dimensional. Average number of iteration steps 
show in different groups AGQISFLA and AGSFLA, as shown in Figure 2-4. 



1 oo 




The number of subgroups 



Figure 2 when G = 100, the average value of the iteration number of different groups of AGSFLA 




5 10 20 25 

The number of subgroups 



50 100 



Figure 3 when G = 1000, the average value of the iteration number of different groups of AGSFLA 
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100 




The number of subgroups 



Figure 4 when G = 100, the average value of the iteration number of different groups of AGQISFLA 

After adaptive grouping, the percentages of the number iterations for each subgroup in total iterations 
are shown in Table 4. 



Table 4 the percentages of the number iterations for each subgroup in total iterations (%) 



Algorithm 


Iterations 


Number of Subgroups 


2 


4 


5 


10 


20 


25 


50 


100 


AGS FLA 


100 


0.14 


1.22 


2.04 


3.32 


5.08 


5.64 


5.78 


76.78 


AGS FLA 


1000 


0.02 


0.08 


0.17 


0.35 


1.87 


3.18 


6.66 


87.67 


AGQISFLA 


100 


0.04 


2.22 


3.62 


2.86 


4.10 


5.18 


6.52 


75.46 



The experimental results demonstrate the AGQISFLA long run time, high capacity optimization 
features, we give the following analysis. 

First, as previously described, the more subgroups, the more number of frogs need to update for each iteration, 
and so the Longer the time of single iteration. Secondly, the individual coding method based on multi-bit 
probability amplitude, if the number of qubits encoded set is n, each iteration subgroup evolutionary times for 
LS, then each iteration, each subgroup of the worst frog are updated LS x n times, which in the extended run 
time, but also greatly improve the number updates of worst frog; at the same time, this by individually adjusting 
the phase of qubits to cycle update individual approach, making individuals more elaborate update, which also 
enhances the solution space of ergodic. Third, in the update strategy based on multi-bit revolving door, and 
draws on the thinking of particle swarm optimization, taking the leading role of subgroups optimal frogs and 
global optimum frog, To some extent, to avoid the tendency to fall into premature convergence; at the same 
time, the use of multi-bit quantum revolving door, one operation can be achieved for all the updates on the 
individual probability amplitude, the characteristics of quantum rotation gates guaranteed probability amplitude 
of the "length" unchanged, effectively avoiding iterative sequence divergence, and thus improve the 
convergence capability. In summary, the adaptive grouping and multi-bit encoding probability amplitude of 
these two mechanisms, at the expense of time in exchange for the ability to optimize, which is consistent with 
the theorem no free lunch. 

IX. CONCLUSION 

In this paper, a quantum inspired shuffled frog leaping algorithm algorithms is presented which 
encoded by adaptive grouping and multi-bit probability amplitude. Frog group individual coding approach is to 
use multi-qubits system in the ground state of the probability amplitude, frog group of individuals update 
method is a multi-bit quantum revolving door. Function extreme optimization results show that under the same 
running time, the optimization ability of proposed algorithm has greatly superior to the conventional leapfrog 
algorithm, thus revealing, adaptive grouping and multi-bit encoding probability amplitude of these two 
mechanisms is indeed an effective way to greatly improve traditional leapfrog algorithm optimization 
capabilities. 
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